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Abstract 

The nucleation rate derived in the classical theory contains at least one 

undetermined parameter, which may be expressed in terms of the Langer 

first-principles theory. But the uncertainties in the accounting for fluctuation 

modes, which are either absorbed into the free energy of a critical cluster or 

not, result in different evaluations of the statistical prefactor and nucleation 

rate. We get the scaling approximations of the nucleation rate for the vapour 

condensation both near and out of the critical range. The results obtained 

deserve the experimental verification to resolve the theoretical uncertainty. 
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The homogeneous nucleation occurred in metastable states has been area of much the- 
oretical and experimental activity since time of Maxwell and Van der Waals. The classical 
nucleation theory was worked out about 60 years ago mainly by Volmer, Becker and Doring, 
Zeldovich, and Frenkel f^-^. Under the treatment of relaxation of the initially homoge- 
neous metastable state as a kinetic process, triggered by the creation of spherical clusters of 
the new, stable, phase as thermal fluctuations in the metastable phase, the classical theory 
provides us the fundamental nucleation rate 

where /q is a pre-exponential factor, AFc is the minimum work needed to create critically 
large cluster of radius R^-, is Boltzmann constant, and T is the temperature of the system. 
The nucleation rate yields the number of viable clusters of the new phase per unit volume 
passing through the critical region per unit time via the equilibrium number of critical 
clusters /o(-Rc) given by the theory of thermal fluctuations. Usually Eq. (0) is rewritten 
in the form 

I = B{R,)UR,)Z , (2) 

containing the size diffusion coefficient B{Rc) and the Zeldovich factor 

7 f d'AFiR) \ 
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This is the principal result of the classical theory. According to the thermodynamic theory 
of fluctuations the multiplicities of clusters of various sizes obeys the Boltzmann (or rather 
Gibbs) distribution 



fo{R) = Co exp 



and the coefficient Co is not determined in the classical theory. To find out the diffusion 
coefficient -B(-Rc) oc 1/ (R — Rc) (d-R/dt)/j=/j^ one has to solve the macroscopic equations con- 
sidering, for instance, the growth of a cluster due to the diffusion flux through its interface. 



The coarse-grained field theoretical approach to the nucleation phenomenon was devel- 
oped by Langer who generalized the earlier works of Landauer and Swanson [Q, and 
Cahn and Hilliard . The prefactor Jq in the theory of Langer is shown to be a product of 
the dynamical prefactor k and the statistical prefactor Qq 

lo = , (5) 

and both prefactors are determined explicitly in the theory. The dynamical prefactor is 
related to the growth rate of the critical cluster, and the statistical prefactor is a measure of 
the phase space volume available for the nucleation. Since the beginning of 70 's, the Langer 
theory has been applied to describe first order phase transition in various systems including 
the vapor condensation [|l^,|ll|], nucleation in binary fluids |]12[, solidification of the melt [jl3 



and hadronization of the quark— gluon plasma [|I4|,|T5|, produced in heavy— ion collisions. 
As was shown by the authors recently |16|, nucleation rate (^) of the field nucleation 



theory may be obtained under the certain assumptions in a form identical to that of the 
classical theory. It means, particularly, that one is able to determine the unknown constant 
Cq in terms of the modern theory. However, in the present paper we show that the difficulties 
in the evaluation of Qq give rise to the uncertainty in the determination of the nucleation 
rate, which can be clearly verified experimentally. 

Our arguments are the following. According to the nucleation theory, both classical and 
modern, the relaxation of the homogeneous metastable state proceeds via the formation of 
clusters of the new, stable, state. The free energy density functional F{{ri}) of the system 
has two local minima in the configuration phase space {rj}, described by a set of collective 
coordinates rji, corresponding to the states of metastable and stable equilibrium. These wells 
in the phase space are separated by the potential barrier related to the formation of cluster 
(or fiuctuation) of critical size. The nucleation rate is determined by the finite probability 
fiux from the metastable configuration to the stable one across the vicinity of point with the 
lowest energy on the potential barrier, the so-called saddle point. 

In the harmonic approximation for F the statistical prefactor appeared in Langer theory 



reads [|T7 

-1/2 



nn = v 
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(6) 
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Here is the volume of the system, Mq and M denote the matrix Mij = d'^AF/drjidrij, 
evaluated at the metastable minimum {770} and at the saddle point, and Ai is the only nega- 
tive eigenvalue of M, associated with the instability of the critical cluster against expansion 
or shrinking. The prime in Eq. indicates that the eigenvalue Ai is omitted. Then, 
the translational symmetry is broken because of the presence of critical cluster. These three 
translational modes with zero eigenvalues give rise to the factor P] proportional to (|Ai|)~'^''^, 
and Eq. (^) is transformed into 

«-m"'n.m"'n(&)"'' 
-"(«)"■• 

From here the uncertainties in the determination of the Qq begin. If the products over /3 
and a are absorbed ^ into the free energy of the metastable state and the saddle point, 
respectively, then Eq. (|) is reduced to 

,,1, _ „ 327r2(tBT)"/2 



If, on the other hand, one will pair [|lOl the eigenvalues A^^-* and A?^ so, that the correc- 



'13 

tions to the free energy difference remain of order R^/V in the limit V —>■ 00, there will be 
four unpaired eigenvalues A^^'^'s at the bottom of the spectrum of eigenvalues 

ii^n(Af)"^(^g)\ (10, 

where the second derivative of the Helmholtz free energy / with respect to the local density Ug 
is inversely proportional to the correlation length ^ squared, d'^f/dn'^ oc The statistical 
prefactor in the latter case turns out to be 
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To make a comparison between the nucleation rates calculated with ^I'q^ and fig^'' one 
should know also the dynamical prefactor k and the minimum work needed to form the 
critical cluster. For simplicity, let us consider the process of vapour condensation at the 
condensation point T, corresponding to a certain pressure pq. If the vapour is overcompressed 
with respect to po, the critical radius is determined by the expression 

^^=A^ = ^' (12) 
A/iAn A/i ^ ^ 

where vi is the volume per particle of the liquid phase, An is the difference in local number 
densities, and A/i is the difference in chemical potentials of the vapour and liquid phase 

A/i = /i„ - /i, = ksT In — . (13) 

Po 

Ratio p/po indicates the degree of supersaturation in the system. We see that in the capillary 
approximation the critical radius scales near the condensation point as (Ine)^^. Since for 
the only negative eigenvalue Ai one gets Ai = —2\\kBT [|16|, where the similarity number 



T8| is given by \z = {4:TTa / ksTy^'^Rc, the statistical prefactor qI^'^ should increase with the 



rise of e, 

nS'^oc(lne)^ (14) 
while f2o drops inversely proportional to the same rate 

n^o^ oc (Ine)-^ . (15) 



For the dynamical prefactor we have ||TT],|T2| 



^ = -^^{lner. (16) 

Here A is the heat conductivity and / is the latent heat. The excess of the free energy due 
to formation of the critical cluster is simply 



AF, A| bo 



keT 3 (lne)2 



(17) 



containing the parameter bo = 16TT{(j/kBT)^vf/3. 

Combining our results we get that the nucleation rate is proportional either to the factor 
ai(e) or to the factor 02 (e), correspondingly, 

= ^ exp (-^) oc (Ine)^ exp ( --^] = a,{e) , (18) 



27r ^ "^V kBTj ' ' (lne)2/ 

- i ^ i-B) - (- (1^) - ■ ^''^ 

Figure |l] depicts the evolution of both ai and 0:2 with the changing of the degree of 
supersaturation e for different values of the parameter Bq. The rapid falloff of the exponential 
at e — » 1 dominates the increase of the statistical prefactor i7g , and although 

limQ^o'^ = and limfij^^ = 00 , (20) 

the nucleation rate in both cases drops to zero near the condensation point and rises up when 
the vapour is over compressed. But the slopes of d//de are different, and this disagreement in 
theoretical predictions may be verified experimentally. Note that in the region near the value 
e = e = 2.7318. . . the slopes are similar. Therefore, the measurements should be carried 
out either with relatively weak supersaturations, or with the strong ones, say e ~ 5, which 
is still below the limit of supersaturation |]TP[ obtained in the experiments on homogeneous 



nucleation in vapour. 

Near the critical point we should take into account the power-law approximations for the 
thermodynamic quantities appeared in the nucleation rates, and the critical indices come 
into play. These quantities are listed in Table | as functions of the standard dimensionless 
variable 9 = 1 — T/T^ for the vapour side of the coexistence curve at T < Tc. 

From Table |I| it follows that the difference of the pressures Ap inside and outside the 
droplet is Ap = AnAfi oc and, consequently, Rc cx 0'^'^~l^'"i and \z oc 0^'^~l^~"i , 

It is worth noting that ^ as 6^/^ at T Tc. Therefore the formulae (^ and ( pH]) 
cannot be applied to calculate the nucleation rate because the criterion ksT /\Xi\ ^ 1 
of the steepest descent method is not fulfilled. To estimate the difference between the two 



approaches, however, one may apply the droplet model of Fisher instead of performing 



of numerical evaluation of the integrals. In the Fisher model the only negative eigenvalue 



Ai is m 



Ai = kBTi9T + 2A|) , (21) 

where r ~ 2.2 is the critical exponent. Although the critical radius determined in the 
Fisher model does not coincide with that of the capillary approach, it has the same scaling 
behaviour, oc 0'^^~'^~^ , in the proximity of the critical point. The free energy of a critical 
droplet is given by 

|i£..(l_3.„^)-|, (22, 

where r^m is the radius of a single molecule. This yields for the exponential factor sXT ^ Tc 

exp (^1^) oc ^3^(^+/^-2^) . (23) 
Substituting Eq. ( pT]) in Eq. (y) we find that near the critical point 

fij,^^ oc , (24) 

and, consequently, 

/(^) = COnsti X ^27+3r(7+/3-2i.) ^ ^^^^^^ ^ ^3.6 _ (25) 

Similar estimations of the statistical prefactor VLq lead to the following parametrizations 

Q!-^^ oc 015-47-4/3 ^ (26) 

and 

= const2 X 0l2--27-4/3+3r(7+/3-2.) ^ ^^^^^^ ^ ^5.1 (27) 

We see again that the resulting power-law approximation of the nucleation rate J^^-* as a 
function of 6 deviates clearly from that obtained for 

In conclusion, we study the effect of absorption of the fluctuation corrections to the excess 
free energy of a cluster into the cluster activation energy. To compare different predictions 



made by the Langer coarse-grained field theory of nucleation for the statistical prefactor Qq 
we choose the liquid— gas phase transition near and out of the critical point. It is shown 
that although the statistical prefactor drops to zero in the first case and goes to infinity 
in the other one near the condensation point far from the critical region, the qualitative 
behaviours of the nucleation rates are similar. However, the power-law parametrizations 
of the nucleation rate as a function of supersaturation are different. The disagreement 
between the predictions is especially noticeable in the region of very weak or, in contrast, 
strong supersaturations. 

The scaling behaviour of the statistical prefactor and total nucleation rate is studied 
also near the critical point in the Fisher droplet model. The qualitative dependences of 
and / on the scahng parameter 9 = 1 — T/Tc are similar in both approaches, while the 
analytic parametrizations are different again: 9^'^ vs 9^'^. The experiments on homogeneous 
condensation of a vapour may certainly resolve this theoretical ambiguity. 
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FIGURES 




TABLES 



TABLE I. Power-law approximation of thermodynamic quantities near the critical point 



Quantity 


Power-law approx. 


critical 




at T ^ Tc 


exponent 


a 




u = 0.625 


An 


(Ano)^^ 


p = 0.5 


A 




7 = 1.0 


I 












A/x 


cx 0^ 
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